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We study different techniques that allow us to gain complete knowledge about 
an unknown quantum state, e.g. to perform full tomography of this state. We fo- 
K^ . cus on two apparently simple cases, full tomography of one and two qubit systems. 

We analyze and compare those techniques according to two figures of merit. Our 
first criterion is the minimisation of the redundancy of the data acquired during 
the tomographic process. In the case of two-qubits tomography, we also analyze this 



O 

^O • process from the point of view of factoris ability, so to say we analyze the possibility 



to realise the tomographic process through local operations and classical communi- 
cations between local observers. This brings us naturally to study the possibility to 
factorize the (discrete) Wigner distribution of a composite system into the product 

iPn ' of local Wigner distributions. The discrete Heisenberg- Weyl group is an essential 

ingredient of our approach. Possible extensions of our results to higher dimensions 

O^. are discussed in the last section and in the conclusions. 



c^ 1 Introduction 



The estimation of an unknown state is one of the important problems in quantum 
information and quantum computational |2]. Traditionnally, the estimation of the 
d"^ — 1 parameters that characterize the density matrix of a single qu(iit consists of 
realising d+1 independent von Neumann measurements on the system. For instance, 
when the system is a spin 1/2 particle, three successive Stern-Gerlach measurements 



performed along orthogonal directions make it possible to infer the values of the 3 
Bloch parameters p^,., Py, and Pz defined by 

{(^x) =Px =ism9 cos ip 

{(jy) = py = 'J sin 9 sin (p (1) 

{cTz) =Pz = lcos9 

Once we know the value of these parameters, we are able to determine unambigu- 
ously the value of the density matrix, making use of the identity 

P(7, 0,ip) = -{I + p^a^ + PyCTy + pzaz) = 2 (^ + ^•^) (2) 

When the qubit system is not a spin 1/2 particle but consists of the polarisation of 
a photon, a similar result can be achieved by measuring its degree of polarisation 
in three independent polarisation bases, for instance with polarising beamsplitters, 
which leads to the Stokes representation of the state of polarisation of the (equally 
prepared) photons. 

Tomography through von Neumann measurements presents an inherent drew- 
back: in order to estimate the (P — 1 independent parameters of the density matrix, 
d + 1 measurements must be realised which means that d"^ + d histograms of the 
counting rate are established, one of them being sacrificed after each of the d+1 mea- 
surements in order to normalize the corresponding probability distribution. From 
this point of view, the number of counting rates is higher than the number of param- 
eters that characterize the density matrix, which is a form of redundancy, inherent 
to the tomography through von Neumann measurements. 

Besides, it is known that a more general class of measurements exists than the 
von Neumann measurements. This class is represented by the Positive-Operator- 
Valued Measure (POVM) measurements 0, of which a reduced subset, the Projection- 
Valued Measure (PVM) measurements corresponds to the von Neumann measure- 
ments. The most general POVM can be achieved by coupling the system A to an 
ancilla or assistant B and performing a von Neumann measurement on the full sys- 
tem. When both the system and its assistant are qurfit systems, the full system 
belongs to a d^ dimensional Hilbert space which makes it possible to measure d"^ 
probabilities during a von Neumann measurement performed on the full system. As 
always, one of the counting rates must be sacrificed in order to normalise the prob- 
ability distribution so that we are left with d"^ — 1 parameters. When the coupling 
to the assistant and the von Neumann measurement are well-chosen, we are able in 
principle to infer the value of the density matrix of the initial qudit system from the 



knowledge of those (P — 1 parameters, in which case the POVM is said to be Infor- 
mationnally Complete (IC). Obviously, this approach is optimal in the sense that it 
minimizes the number of counting rates (thus of independent detection processes) 
that must be realised during the tomographic process. 

As it was shown in j4j, the PVM approach to tomography can further be opti- 
mised regarding redundancy. Optimality according to this particular figure of merit 
is achieved when the d + 1 bases in which the PVM measurements are performed 
are "maximally independent" or "minimally overlapping" so to say when they are 
mutually unbiased (two orthonormal bases oi a d dimensional Hilbert space are said 
to be mutually unbiased bases (MUB's) if whenever we choose one state in the first 
basis, and a second state in the second basis, the modulus squared of their in-product 
is equal to 1/d). It is well-known that, when the dimension of the Hilbert space is 
a prime power, there exists a set of d + 1 mutually unbiased bases [U El IH]. This 
is the case for instance with the bases that diagonalize the generalised Pauli oper- 
ators [HIIZI. Those unitary operators form a group which is a discrete counterpart 
of the Heisenberg-Weyl group, the group of displacement operators [Hj, that present 
numerous applications in quantum optics and in signal theory p. 

A discrete version of the Heisenberg-Weyl group [TU] also plays an essential 
role [lij in the derivation of so-called covariant symmetric-informationally-complete 
(SIC) POVM's. Such POVM's are intimately associated to a set of d^ minimally 
overlapping projectors onto pure qudit states (the modulus squared of their in- 
product is now equal to 1/^/d + 1). 

We shall compare the respective merits of PVM and POVM tomographic pro- 
cesses in the cases of one and two qubit systems and focus on the factorisability of 
the latter. This question leads us to study the factorisability of discrete Wigner qua- 
sidistributions which appear to be a very natural tool in the context of one qubit 
and (factorisable) two qubit tomography. Although in the continuous case, the 
factorisation-property is de facto fulfilled, in discrete dimensions this is not a trivial 
question at all. This why this question raised recently a lot of interest (Refs. ^2 ^o 

I23IISll2ni)- 

Generalisations to higher dimensions are discussed in the last section and in the 
conclusions. 



2 Tomography of a (single) qubit system. 



2.1 Optimal PVM approach. 

The aforementioned traditional approach to tomography for a spin 1/2 particle, that 
consists of three successive Stern-Gerlach measurements performed along orthogonal 
directions X, Y and Z is optimal among PVM tomographic processes because the 
3 corresponding bases are MUB's [SHHEnj. Considered so, the traditional approach 
to spin/polarisation tomography is optimal. Actually, the a operators plus the 
identity constitute a discrete counterpart of the displacement operators. Formally, 
they can be defined as follows: cxij = w (—)*■' X]fc=o(^)'^'"'l^ + ^("^od. 2)) (/c|, where 
the labels i,j,k can take values or 1. It is easy to check that, up to a global 
sign that we presently keep undetermined, (To,o = Id.,aifi = <Jx-,<y\,\ = cry and 
<7o,i = <^z- This set is orthonormal regarding the Trace-norm: ^Tf-f^lf^j = Sij (here 
d = 2) and, like the displacement operators in the continuous case, its elements 
constitute a complete basis of the set of linear operators of the Hilbert space on 
which they act. The Bloch parameters are seen to be in one-to-one correspondence 
with the qubit Weyl distribution (defined by Wij = {l/2)Tr.{p.aij)), which consists, 
in analogy with its continuous counterpart ^^ , of the amplitudes of the development 
of the density matrix in terms of the (qubit) displacement operators [9 : wofi = 
l/2,wifi = Px/2,wi^i = Py/2, wq^i = Pz/2. These properties can be generalised 
to higher dimensions |7j, and, when the dimension is a prime power, each PVM 
measurement in one of the d + 1 MUB's leads to the estimation of a set of (i — 1 
parameters of the Weyl distribution. The measurements performed in different 
MUB's are independent and the d+1 subsets oi d— 1 amplitudes obtained so form 
a partition of (the set oi d^ — 1 independent parameters of) the Weyl distribution 
and provide a complete tomographic information about the unknown qut/it state of 
the system. 



2.2 Optimal POVM approach. 

It has been shown in the past, on the basis of different theoretical arguments fH 
EH 123 5 that the optimal IC POVM is symmetric in the sense that it is in one- 
to-one correspondence with a tetrahedron on the Bloch sphere. Intuitively, such 
tetrahedrons homogenize and minimize the informational overlap or redundancy 
between the four histograms collected during the POVM measurement. Some of 



such tetrahedrons can be shown to be invariant under the action of the Heisenberg- 
Weyl group which corresponds to so-called Covariant Symmetric Informationnally 
Complete (SIC) POVM's HIJ. 



Let us now briefly describe how such a POVM measurement could be realised 
experimentally (we were able recently J2Hj to implement this POVM measurement on 
a two qubit NMR quantum computer [21! )• Let us suppose that we wish to estimate 
the three parameters 7, ^,and y? necessary in order to describe the unknown state of 
the qubit a. An ancilla is added to this device as qubit b to form a extending system. 
This device is initially prepared in the state: pj„ = Pa® |0) (0|^. This state differs 



according to different input qubits a. In virtue of the Stinespring-Kraus theorem [2? 
quantum operations are related to the unitary transformations, a property that we 
shall now exploit by letting the entire system evolve under unitary evolution U 
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where a = Jl + l/-\/3, /? 



1/v^. 



By measuring the full system in a basis that consists of the product of the a and 
h qubit computational bases, we are able in principle to estimate the four parameters 
enlisted on the diagonal of the following matrix: 
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such a POVM measurement is informationally complete due to the fact that the 
PqO) -Poi) -PiO) -Pii are in one-to-one correspondence with the Bloch parameters p^,, Py, 
and p^ as shows the identity 
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Actually, Pqo is the average value of the operator (|)(cro,o + (-7^)(o"i,o + <^o,i + 
(Ti^i)) which is the projector onto the pure state |0)(0| with |0) = a\0) + /3*|1) 
and a = Jl + 4=, /?* = e^^l — 4=. Under the action of the Pauli group it 
transforms into a projector onto one of the four pure states aij\(j)); i,j : 0,1: 
^^M{<f>W^,J = (i)((l - 73)^0.0 + ipi^ii=oi-y'-'-'crk,i)) The Signs i-y-^-^-'^ re- 
fleet the (anti) commutation properties of the Pauli group. So, the four parame- 
ters Pij are the average values of projectors onto four pure states that are "Pauli 
displaced" of each other. The in-product between them is equal, in modulus, to 
l/y/S = l/y/d + 1, with d = 2. This shows that this POVM is symmetric in the 
sense that it is in one-to-one correspondence with a tetrahedron on the Bloch sphere; 
as this tetrahedron is invariant under the action of the Heisenberg-Weyl group it is a 
Covariant Symmetric Informationnally Complete (SIC) POVM jTT]. One can show 
im 1211 I2S] that such tetrahedrons minimize the informational redundancy between 
the four collected histograms due to the fact that their angular opening is maximal. 

It is worth noting that this POVM possesses another very appealing property 
which is also true in the qutrit case but not in dimensions strictly higher than 3 
|29j : the qubit Covariant SIC POVM is a direct realisation (up to an additive and 
a global normalisation constants) of the qubit Wigner distribution of the unknown 
qubit a. Indeed, this distribution W is the symplectic Fourier transform of the 
Weyl distribution w (already defined by the relation Wij = (l/2)Tr.(p.crj j)) which 
is, in the qubit case, equivalent (up to a relabelling of the indices) to its double 
qubit-Hadamard or double qubit-Fourier transform: 



W,,i = (1/2) Y: (-)^-'-^'-V, = ((l/v^) E(-)^')((l/V2) Ei-r'->^,r (6) 

i,j=0 i=0 i=0 



This expression, originally derived by Wootters in a somewhat different form 
[Oj . is a special case of an expression for a Wigner distribution derived by us in 
prime power dimensions jHOl- One can check that Pk,i = {'^/V^)Wk,i + (1 — l/"\/3)/4. 
The symplectic Fourier transform is invertible so that once we know the Wigner 
distribution, we can directly infer the Weyl distribution or, equivalently, the Bloch 
vector of any a priori unknown quantum state. It is worth noting that in order 
to measure the coefficients Pk,i{Wk^i) we must carry out a measurement on the full 
system (the unknown qubit plus the ancilla), which is the essence and novelty of 
entanglement-assisted quantum tomography 



It has been shown that the discrete qubit Wigner distribution directly gener- 
ahses its continuous counterpart [3T] in the sense that it provides information about 
the locahsation of the qubit system in a discrete 2 times 2 phase space [12]. For 
instance the Wigner distribution of the first state of the computational basis (spin 
up along Z) is equal to VFfc,/(|0)) = {l/2)6kfi, which corresponds to a state located in 
the "position" spin up (along Z), and homogeneously spread in "impulsion" (in spin 
along X), in accordance with uncertainty relations (see Ref. [20] for an enlightening 
discussion of discrete uncertainty relations in connection with the Wigner distribu- 
tion). Similarly, the Wigner distribution of the first state of the complementary 
basis (spin up along X) is equal to Wk,iiil/\^){\0) + |1))) = (l/2)5/,o- 

We arrived to our formulation of the Wigner distribution 30 by deriving a 
solution of the Mean King's problem j331 1^ EH] , which is not astonishing. Indeed, 
this problem consists of ascertaining the value of the spin of a spinor prepared at 
random in the X, Y and Z bases. The connection between the Wigner distribution 
and the Mean King problem is the following. Its solution consists of entangling 
the qubit to another qubit and to measure the full system in a well-chosen quartit 
basis in such a way that each detector fires with a probability equal to the Wigner 
distribution of the first qubit. Therefore, knowing to which basis the initial states 
belong and knowing which of the four "Wigner" detectors would fire, we are able to 
infer what is the value of their spin component [3211201. For instance, when the spin 
is prepared in the Z basis and that the detector corresponding to Wi^i {i = 0, 1) 
fires, we could infer that the initial spin was the state |1) (spin down along Z). 



3 Tomography of a two qubit system. 



3.1 Optimal PVM approach. 

The results relative to the qubit case can be generalised to multi-qurfit systems 
whenever d is a prime power [3] , so to say to the case d = p^ with p prime and m 
a positive integer. In particular, in the present case {p = m = 2), we can define 
generalised displacement operators according to the definition given in [7j (see also 
Ref.[2ni): 



yj ^ J24'''''''''''''^\k®G^){k\■,^,J : 0...d-l, (7) 

fe=0 



where 7g is the pth root of unity: jg = e*-^'^/^, the Galois addition Q)g and the 
Galois multiplication ©^ are defined by the tables given in appendix. The Ga- 
lois addition is by definition equivalent with the addition modulo p componen- 
twise. Concretely, this means that if we write the labels (0,1,2,3) in a binary 
form (0=(0,0),1=(0,1),2=(1,0),3=(1,1)), the Galois addition is defined as follows: 
if (i = {ia,ib) and j = {ja,jb)), then i ®g j = {ia ®mod2 ja,ib ®mod2 jb)- The 
Galois multiplication is distributive relatively to the Galois addition, moreover it 
is commutative, and there is no divider of (the neutral element for the addi- 
tion) excepted itself. The algebraic structure that is defined by these require- 
ments is a field; in the present case, the multiplication table is uniquely defined by 
these requirements, and by the definition of the addition: ©g x = x ©g = 0; 
1 ©G x = X ©G 1 = a;;2 ©G 3 = 3 ©G 2 = 1 and 2 ©G 2 = 3. 

We can write the 16 4 times 4 displacements operators defined by Eqn.(|7j) as the 
identity plus 5 sets of 3 operators that are defined as follows. The first set consists of 
the operators V^' with z = and j = 1, 2, 3, while the other sets consist of the opera- 
tors Vi with i = 0, 1, 2, 3, and correspond to the respective choices / = 1, 2, 3, 4. 
By direct computation, one can check that we obtain so the five sets {a", a^, a'^.a^}, 
{<, al <.a^}, {c^^^ cr^, cr^^a^}, {<.a^, a«.o^,^ <.a^} and K.a,^ a^^.a^ a^.a^} (up to 
irrelevant global phases). The operators that belong to each of these sets commute 
with each other; moreover, it is possible to multiply each of them by a well-chosen 
phase in such a way that the 3 operators of each set form (together with the identity 
operator) a commutative group. The bases that simultaneously diagonalize all the 
operators of such sets are unambiguously defined and are mutually unbiased jHl El 
relatively to each other. All the properties that we described in the qubit case are 
still valid in the present case: by performing von Neumann measurements in those 
5 MUB's it is possible to estimate 15 parameters (5 times (4-1) probabilities) that 
are in one-to-one correspondence with the Weyl distribution, or equivalently with 
the coefficients of the density matrix of the system. 

The problem is that only the three first bases are factorisable (the common 
eigenbases of the two last operators are actually maximally entangled [HI 13 IH7]). 
This is not astonishing because by taking products of mutually unbiased qubit bases 
we can at most construct 3 MUB's. Of course, it is still possible to obtain full 
tomographic information about the system by measuring each qubit component 
in the three MUB's (along X,Y and Z), due to the fact that the displacement 
operators are factorisable into a product of local displacement operators (a very 
general property, also valid in dimension p™ with p an arbitrary prime and m an 
arbitrary positive integer [7,). The problem is that this requires to establish 3^.2^ = 
36 count rates in order to estimate 15(-|-1) parameters, a situation which is far 

8 



from being optimal. To conclude, it is clearly impossible to conciliate optimality 
and factorisability in the two qubit case (as well as in the two qudit case, when 
d = [p"^)^) because factorisable MUB's are necessarily products of local MUB's 
{vcP = yd.yd) and, by taking products of mutually unbiased bases of the p^ 
dimensional Hilbert spaces associated to the two components of a bipartite system 
of dimension d = (p™")^, we can at most construct p"^ + 1 (factorisable) MUB's 
which, for all prime power dimensions, is strictly smaller than d + 1 = {p^Y + 1- 



3.2 Optimal POVM approach. 

As in the two qubit case, it is easy to show on the basis of a simple argument 
that it is impossible to conciliate optimality and factorisability of the POVM or 
entanglement-assisted tomography in the two qubit case (as well as in the two qufiit 
case, when d = (p"^)^). The reason therefore is that when d = di.d2, and that a 
POVM can be splitted into a product of two SIC POVM's which means that we 
couple the di (^2) dimensional subsytem to a di (^2) dimensional assistant we find 
d"^ = d\.d\ projectors onto factorisable pure states of which the in-products are most 
often equal in modulus to [1/ ^/dl^^).{l/ y/d^^^) but not always; sometimes this 
in-product is equal to {1/yJdi + 1) or (l/y32~+T), so that what we get is not a SIC 
POVM, but only an IC POVM which may not be optimal. 

Nevertheless, the number of counting rates necessary in order to realize a to- 
mographic process by a factorisable POVM measurement is optimal and equal to 
d^ = d\.d\, so that this technique can be considered as a good compromise: not 
fully optimal from the point of view of redundancy but at least factorisable which is 
very appealing regarding the experimental realisability of the tomographic process 
in the case of separated subsytems (for instance in the case of two separated photons 
entangled in polarisation, a common situation in the laboratory). In the case of two 
qubits, a factorisable SIC POVM is nearly equivalent (actually, it is equivalent up 
to a simple renormalisation of the counting rates) to a direct measurement of the 
average values of the 16 products of the local (qubit) Wigner operators defined in a 
previous section (Eqn.®. As we are interested in the question of the factorisability 
of the tomographic process, it is very natural in the present context to ask the fol- 
lowing question: Is the product of the two Wigner distributions of the subsystems 
of a bipartite system equivalent to the Wigner distribution of the full system? We 
shall provide certain answers to this question in the following section. 



4 Factorisability of the discrete Wigner distribu- 
tion. 

4.1 Candidates for the Wigner distribution. 

Before we discuss their factorisability, it is necessary to provide an operational defi- 
nition of the discrete Wigner distribution for a finite-state system. We derived such 
a definition in R.ef . [HH] where we showed that the recipe for constructing a Wigner 
distribution associated to a qudit system, with d = p"^, was the following: 

i) Let us split the set of cP d times d displacement operators defined by Eqn.(|Zj) 
into the identity plus d + 1 sets oi d — 1 operators that are defined as follows. The 
first set consists of the operators V/ with i = and j = 1, ...,d — 1, while the other 
sets consist of the operators Vi with / = 1, ...,d — 1 and correspond to the 

respective choices I = 1, ...,d. In the qubit case, each set corresponds to one of the 
Pauli operators. In the two qubit case, the list of these sets was explicitly given 
in the section 13.11 It can been shown |SH] that all the operators from a same set 
commute with each other. 

ii) Let us multiply each operator of a set by a well-chosen phase in such a way 
that the set of "renormalised" operators together with the identity forms a finite 
group (with d elements). 

It is shown in Ref. [7] that there are, for each of the d+1 sets, d possible choices of 
phases that satisfy this constraint. Moreover these choices are equivalent to simple 
relabellings (in fact translations or Galois-additive shifts) of the indices of the states 
of the MUB in which the operators of the group are simultaneously diagonal. 

In the same reference two possible choices are explicitly given, that correspond 
to the odd and even dimensional cases: 

In the odd case, we showed that one among the possible phase choices led to 
the following definition of the renormalised displacements operators associated to 
the ith group (denoted U}, with i : 1, ...,d, I : 0, ...,d — 1): 

where /g represents the inverse of the field (Galois) multiplication and 2 = 1 (Bg 1- 
Actually, this choice is particularly attractive and elegant for several reasons and is 

10 



uniquely defined once we know the operation tables of the field with p"^ elements. 
We shall show that it helps us to answer positively to the problem of the factorisation 
of the Wigner function for bipartite systems in odd prime power dimensions in a 
next section. 

In the even case {d = 2™) the explicit expressions for the phases are less easy to 
manipulate, essentially due to the fact that 1 (Bo 1 = and that we may not divide 
by 0. Once again, there are p"^ (2™ in this case) possible ways to determine the 
phases f/f/VJ which are equivalent, up to a relabelling of the states of the 

corresponding MUB. 

A possible choice for the phases was shown to be 

where the coefficients /„ are unambiguously defined by the p-ary (here binary) ex- 
pansion of /, / = Y.T=o ^n2", while n' is the smallest integer strictly larger than n 
such that /„/ 7^ 0, if it exists, otherwise. 

For i = 0, a possible choice of phases corresponds to the relation Uj^ = V^, in 
the even and odd cases. 

It is worth noting that in both cases the phases are square roots of integer powers 
of gamma, 

(f//7v;^^"'^®^')2 = ^|G({i-i)0G/0GO_ As a consequence, it is easy to show (see 
appendix) that (f//)-^ = (f//V = U^^i 

iii) In Ref . [^ . the cP Wigner operators are defined as follows: 



"• m,n=0 

Introducing the more convenient notation Um,n = (7g ' ^"^ J^^m i^^^ Um,n operators 
are equivalent with the Uj operators previously defined, up to a mere relabelling), 
we get 

Win,.) = ]i E 7l^^^®°"®°^^®^'"^r.,n). (11) 

"■ m=l,n=0 

11 



It is worth noting that the Wigner (W) operators defined by Eqn. ljlUj) are Her- 
mitian, due to the identity [/^„ = U~\ = UQ^rn,eGn- In appendix, we prove that 
they are "acceptable" candidates for discrete Wigner distributions according to the 
criteria introduced by Wootters in his seminal paper of '87 ^Hj, which means that 
(a) their Trace is equal to 1, (b) they are orthonormalised (to d) operators (under 
the Trace norm) so that they form a basis of the set of linear operators, (c) that if 
we consider any set of parallel lines in the phase space, the average of the Wigner 
operators along one of those lines is equal to a projector onto a pure state, and the 
averages taken along different parallel lines are projectors onto mutually orthogonal 
states. In appendix we prove that the last relation is valid in odd and even prime 
power dimensions as well and we also show that the sets of d orthogonal states that 
are associated to different directions (there are d -\- 1 non parallel directions in the 
d times d phase space) form MUB's (in accordance with the prediction made in 
Ref-dH). 

The Wigner distribution is then nothing else than the set of the d"^ amplitudes 
that we obtain when we develop the density matrix of a qudit (with d = p^) in the 
basis provided by the Wigner operators: 

W(n,i2) = Tr.{p.W(^,,^i,,)) 

As a consequence of the property c the marginals of this distribution along 
any axis of the phase space are equal to transition probabilities to states of the 
corresponding MUB. 



4.2 Factorisability of the two and three qubit Wigner dis- 
tribution. 

In Ref ■ PU] ■ we showed that W^(ji,i2) = (^/i^)~"^W^(o,o)^i^ ? and we also showed in Ref. [7j 
that the displacement operators always factorise into products of local displacement 
operators. Therefore, in order to prove the factorisability of the Wigner operators, 
it is sufficient to prove that 14^(0,0) factorises. As we mentioned in the previous 
section, the Wigner operators are not uniquely defined according to our definition. 
For instance, in the qubit case we are free to change the sign of the operators 
ax, (Ty, and az according to our convenience in the definition (jH)). This change of 
sign is equivalent to a relabelling of the states of the associated MUB's which has 
no important physical consequence. Therefore 2^ acceptable Wigner distributions 
exist in the qubit case. In the two qubit case, we are free to choose arbitrarily the 
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sign of two operators in each of the 5 famihes of 3 operators that we defined in the 
section im which corresponds to 4^ acceptable Wigner distributions. There are thus 
(2^)^ possible products of two qubit Wigner distributions and 4^ possible two-qubit 
(quartit distributions). Once we have chosen the signs of the three sigma operators 
{x,y,z) of the first qubit, the corresponding signs at the level of the second qubit 
must be the same an even number of times (0 or 2 times) in order that the product of 
the Wigner distributions is an acceptable two qubit distribution. For instance, if we 
choose the phase + for the a sigma operators (x, y, z) and the phase — for the h sigma 
operators (x, y, z), we obtain the products {+cr", — cr^, — cr^.a^}, {+cr°, — o"^, — o'^.cr^}, 
{+^^ -^l -^r^'y]^ {-<•<, -^r< -^>l] and {-<.cr^, -a^-al -a^}. Each of 
these families (together with the identity) forms a commuting group, so that the 
corresponding operator 14^(0,0) , which is the sum of those 16 operators factorises: 



1 

'(0,0) 



IVro.o^ = -^{+Id. + a: - o\ - al.o\ + < - a^ - a^.a^ + 



^y-S 






h^^Id.^^ol^ol^al)\ 



-i+Id.^ + < + < + a:).-i+Id.'' ~al-al- a^). (12) 



As a consequence, all Wigner operators factorize too. There are obviously 2^.4 = 
2^ similar ways to derive factorisable two-qubit Wigner distributions, so that fifty 
percent of the products of qubit Wigner distributions provide an acceptable quartit 
Wigner distribution (among the 4^ acceptable Wigner distributions). 

It is now easy to show that a three qubit Wigner distribution never factorises 
into the product of three qubit Wigner distributions. Essentially this is due to the 
fact that it is impossible to find three ordered triplets of plus or minus signs that 
would be the same an even number of times (0 or 2 times) TWO BY TWO. 

It is important to note that these results are still valid in the approach followed 
in Ref. |17j. This (more axiomatic and geometrical) approach is more general be- 
cause it allows more flexibility in the way to attribute one element of the Galois field 
to one label of the basis states of the computational and dual bases. Implicitly in 
our approach, the muple of integers comprised between and p — 1 assigned to each 
computational state expresses the development of the corresponding element of the 
Galois field in a field basis (field bases are defined in Ref.|E|) that contains 1 (the 
neutral element for the multiplication) as first element. Similarly, our choice of the 
dual basis is such that the MUB associated to diagonal straightlines factorises into 
local MUB's associated to local diagonal straight lines [7j. Anyhow, the MUB's are 
exactly the same in both approaches, independently on the choices of field bases, 
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because the splitting of the Heisenberg-Weyl into commuting sub-groups is unam- 
biguously defined up to relabellings once we associate ox operators to horizontal 
translations (shifts of the labels of the first MUB) and oz operators to vertical 
translations (shifts of the labels of the second MUB). In our approach, different 
choices of phase conventions for the subgroup of the Heisenberg-Weyl group lead 
to different labellings of the corresponding MUB while in the approach of Ref. f^I 
this labelling is arbitrarily imposed. Anyhow, there are in both cases (c?)'^"'"^ possible 
phase conventions in dimension d once the field bases are chosen. There are thus 
in both approaches [(^Y ^^ possible Wigner operators in dimension rf^ and (^^'^^^' 
possible factorisable Wigner operators, and those operators are the same because 
they can always be written (up to additive and normalisation constants) as sums of 
projectors onto (i+ 1 states from different MUB's [17j, and the MUB's are the same 
in both approaches. 

One can check for instance that one of the two factorisable Wigner distributions 
derived in Ref.[T7j in the case d = 4 corresponds to choosing the phases (-(-, -(-, -(-) 
for the a sigma operators (x, ?/, z) and the phases (-I-, — , +) for the h sigma operators 
(x, y, z). The second one is obtained similarly but with a permutation of the roles of 
a and h. The 32 operators that we derived can be obtained from those two operators 
by performing at most one local rotation of 180 degrees around the X^Y ox Z axis. 



4.3 Factorisability of the two qurfit Wigner distribution with 
d = p^ and p odd. 

We shall now prove the following result: provided we apply in odd prime power 
dimension the particular phase convention (jH} and define the field with (P = (p"*)^ 
elements as the quadratic extension of the field with d = p^ elements (this tech- 
nique was succesfuUy applied by us in the past in order to solve the Mean King's 
problem in prime power dimensions [201 ), the Wigner distribution of a bipartite p^™- 
dimensional system naturally factorises into a product of local Wigner distributions 
for the two p'"-dimensional subsystems. 

Before we do so, we must define the quadratic extension of a field. Let us denote 
la {ib) the elements of the field with p™ elements associated to the a [b) subsystems. 
Their quadratic extension is a field associated to the bipartite system a — b the 
elements of which we denote {ia ,4). As always, its addition (denoted ©©g) is the 
addition componentwise: 
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{la ,ib) © ®GiJa,jb) = {ia ®G ja h ®G3b)- 

In particular, 2 = (1„ ,0^) © ®g{K. %) = (!« ©g la, 06) = (2„ 0^) 

All what we need to know about the extended multiplication rule (denoted 00g) 
is that it is commutative and distributive relatively to the extended addition, that 

{ia,0)QQG{ja,0) = {laQGja ,0) , 

{ia,0) © ©g(0, J6) = (0,ia QGJb ), and finally that 

(0, ib) © ©g(0, J6) = («b ©G jb ©G R, ib ©G jb ©G Q), with R and Q elements of 
the field with p™ elements, and R different from (otherwise this extension would 
not form a field). Those properties are very similar to those met in the case of the 
complex field which is the quadratic extension of the real (infinite) field. 

It is instructive to note that as 2 = (2(j, 0^), its inverse (relatively to the extended 
multiplication) is equal to (2~^, Ob), where 2~^ represents the inverse of 2 in the non- 
extended field with p^ elements. 

The Wigner operators of the composite system can now be written according to 
Eqn. ()lU|) . and their factorisability is easily established: 



1 d— 1 

Ti/ -•- V^ ^©©G(«i.«J)00G(«a,ni,)eeG(«a:«i)00G(ma,r?i6)/ ({ma,mb)00G(na,n6)xiTr(na,n6) 

^{iil,il)^i^a,^,)) = ^ 1^ ^G 17g )^^{m,,n,,) 



_ 1 Y^ 00G(ii,«J)00GK,ni,)©©G(ia,«b)00G(ma,mi,) {,{ma,rrn,)Q&G{na,rn,)/ / G'i)\ 

- -^ 1^ 1g 1g )■ 

(ma,mf,),(na,n,,)=0 

E 7^"'"''^^®®"^""'"^^^|(A:a,fc.)©©G(m„mb))((fc.,fc,)| 

_ 1 \r^ eG(«i0Gna©G«j0G-R0G"6)©G(«o0GmaeG«i0G^0Gm6) (ma0Gna©G(mb0GR0G»T-6)0G2~^ 

- ^ 2^ 7g 7g 

{ma,mi,,na,rn,)=0 

E 4''^^on.)^-(',eaRean,))^j^^ ^^ m^){K\\h ©G m,)(A:,| 

ka,kf, = 
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5 Conclusions. 



At first sight, the tomography of single and two qubit systems seems to be a trivial 
question. From the previous treatment we see that if we analyze the problem at the 
light of two criteria (optimality in the sense of minimal redundancy and factoris- 
ability) the problem is surprisingly rich. In particular it motivates the interest of 
studying the possibility to factorise the Wigner distribution of a discrete composite 
system, a question that recently attracted an increasing interest f Refs. [T^ to 



Besides the question of tomography of composite systems, the main reason 
therefore is that the phase space structure [121 HE] of composite systems is not 
necessarily factorisable (for instance the two qubit straight line of slope 2 of the 
16 dimensional phase space contains the 4 couples (0, 0)a,b, (1, 2)a,fe, (2, 3)a,b and 
(3,1); it is obviously not the Cartesian product of two qubit straight lines because 
{Oa,Oa) = {0,0)a,ila,2a) = (0, 1)„ (2„, 3j = (1,1), and (3„1J = (1,0), ). Ac- 
tually, the existence of non- factorisable lines is directly related to the existence of 
entangled, non-factorisable MUB's and is an unavoidable feature of composite sys- 
tems of prime power dimensions, for dimensional reasons similar to those that we 
explained at the end of section 13.11 This motivates the quest for global (non nec- 
essarily factorisable) phase space approaches, an example of which is provided by 
our displacement operators and our Wigner distribution. Although the phase space 
of the composite system is not factorisable (it is not the Cartesian product of the 
phase spaces of the composite systems), it could be that Weyl or Wigner operators 
nevertheless factorize, which is the object of the present paper. Of course, the prod- 
uct of local Wigner functions can always provide a full tomographic representation 
of the state of a composite system, but this picture is naturally linked to a Cartesian 
splitting of the full phase space, which is not as rich as a global description of the 
full phase space. The affine structure of the full space phase (which is intimately 
related to the underlying field with p™ elements) is lost whenever we replace it by 
the Cartesian product of the local phase spaces (as is done in Ref.jTHj) as showed 
our example at the beginning of this section. 

There exist several approaches to the problem (Refs.[T2] to j22l HSl dl 120] ) , and 
most often those approaches have a pronounced geometrical flavour in the sense 
that they aim at deriving potential candidates for the Wigner distribution from 
general considerations about the structure of the d times d phase space (an excellent 
survey of the phase space approach is given in the introduction of Ref.^J). Our 
approach is slightly different from the beginning because we postulate from the 
beginning (and this is an educated guess) the "algebraic" expression of the Weyl and 
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Wigner operators (or phase-point operators following the terminology introduced by 
Wootters in [13]). It seems nevertheless that our approach captures the essential 
features of the more general, geometric, approach. For instance, it is also true 
in our approach that straight lines of the phase space correspond to states, and 
that the states associated to parallel lines form orthogonal bases, while different 
orientations correspond to MUB's [inillZI. Our results about the factorisability of 
Wigner distributions are still partial results, and they directly raise a question the 
answer of which is out of the scope of the present paper: 

Is it possible to factorise the Wigner distribution of a (2™)^ dimensional system 
into a product of two (local) Wigner distributions of (2™) dimensional subsystems 
when m is strictly larger than 1? 

Finally, it would also be interesting to investigate the factorisability of Wigner 
operators in odd dimensions (in which case it has been shown that many results valid 
in odd prime power dimensions can be transfered nearly integrally [ ^ I39 | 1 ^ HT] ) . 
For instance the definitions of the Weyl ((Zj) and Wigner operators (|Tn|) are still 
operational when we replace the Galois operations by the modulo d operations and 
the pth root of unity 7g by the dth root of unity. Factorisation is still possible in this 
case. For instance when d = 15 = 3.5, we can write rua^b = 5.ma + 3.mb{modulol5) , 
where < ma,b < 14, < rua < 2 and < rub < 4. Then factorisation is ensured by 
the identities ma,b + na,b{jnodulol'i)) = 5.{ma + na{modulo3))+3.{m,b + nb{modulo5)) 

and exp is =exp s .exp s 

Prime power dimensions remain exceptional anyhow because Galois fields with 
d elements are known to exist only when d is a prime power. Our guess is that a 
set oi d+ 1 MUB's only exists in prime power dimensions jl2], due to the fact that 
they seem to be closely related to the existence of finite fields and finite projective 
spaces O Cni HH] (finite affine spaces with d'^ elements do not exist when d = 6 [Hj 
or 0? = 10 ^Hj and it is conjectured that they exist only when d is a. prime power ), 
but this question is still open. 

Last but not least, it is worth to mention the approach to the problem that was 
developed by Rubin and Pittenger |Tni 1201 • This approach is also algebraical but the 
authors make use of sophisticated techniques like field extensions of arbitrary order 
m (in the treatment of odd dimensions we limited ourselves to quadratic extensions), 
which enabled them to prove that it is possible to find, when the dimension is an 
odd prime power {d= p™ with p prime and ODD), Wigner distributions that would 
be the product of m local Wigner distributions. Our result of the section H7S1 allow 
us to answer positively whenever ?tt, = 2", n = 1, 2, 3.... 
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The authors also estabhshed in Ref. [^ the factorisabihty of the two qubit 
Wigner distribution but their method did not allow them to treat systems composed 
of more than two qubits. As we showed in this paper in the even case [p = 2) we 
know that factorisabihty is not possible when m = 3 which establishes a clear 
distinction between the even and odd (prime power) dimensions. 

Another adavantage of their approach is that they are able to estimate the 
degree of separability of MUB's states jTHj in arbitrary prime power dimensions 
(this reference was kindly drawn to my attention by the authors). It seems that our 
approaches are closely related (for instance, in order to establish the separabihty 
of Wigner distributions [20], the authors also made use of the relative freedom in 
the assignment of phases to commuting operators of a same subgroup). In the 
introduction of Ref.^|, Rubin and Pittenger wrote, relative to the approach of 
Wootters et al. that "Although the motivations of the two approaches appear to 
he quite different, they require the same mathematical tools and appear to lead to 
the same results. An interesting question is the interrelationship between the two 
approaches. " This is certainly true concerning our approach and theirs. 
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Appendix 1: "Acceptability" of the Wigner operators. 

Let us consider the Wigner operators defined by Eqn. (jlO|) . We shall now show that 
(a) their Trace is equal to 1, (b) they are orthogonal with each other and normalised 
to d (under the Trace norm) (c) that if we consider any set of parallel lines in the 
phase space, the average of the Wigner operators along one of those lines is equal to 
a projector onto a pure state, and the averages taken along different parallel lines 
are projectors onto mutually orthogonal states. 

Three identities, that were derived in Ref . [Zj , will be helpful in the establishment 
of the proofs: 

E7r"^^ = A,o (14) 

i=o 

7g ■ 7^ = 7r"^'^ (15) 



\7y^ = 7''^^'^-'^vSf. (16) 

With the help of the identity (jl4|) and on the basis of the definition ((7j), we 
get that tr-Vi = Z]fc~feLo7G ^ ^ ^^^k(BGi,k'^k,k'= d.6ifi.6jfi. It is easy to show that 
the V operators defined in ((Tj) are unitary with (V/)+ = (1^/)"^ = 7®G(^0Gi)ye^G\ 
Making use of the composition law (|16|). it is easy to show that the V operators are 
orthogonal relatively to the trace norm: and that tr.V/ = N.6ifi.6jfi. 

Hence, we can derive the property (a): 
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1 v^d-1 ^,GGilQGn®Gi2QGm( imeGn)\^rp /t^„n 
■^ 2^m,n=0lG \IG J^-'-'-K^m) 

1 X^d-1 .,eGilQGn®Gi2eGm( {meGn)^^ ^ jr r \_lY^d-l ^_1 



In order to prove the property (b), we should firstly note that {Ui)~^ = {UiY = 
Uq^i a direct consequence of the Eqn. (|TT)|l and of the identity (Ui /Vi ''^ y = 
1g° ' '^ • Besides, the U operators, like the V operators are orthogonal 

relatively to the trace norm, so that Tr.{U}^^.Um',n') = d.Sm,m'-Sn,n'- 

Henceforth, Tr.{Wl^^^yW^,^,,)) = ^(^^-i^_^^^^eon0aneo..0a..^eo.i0anea.^0o^^ 
=S{ii,i')S{i2,i'), where we applied twice the identity (fT^ . 

In order to prove the property (c), it is useful to recall the transformation law 
of the U operators that was established in Ref . jSDj : 

( (eG((»-l)0G™©G")0G'")\^ 
Um,n{^) = , eG((i-l)©G'"0G'")\^/ ©g('"Og")>^ ©Gfi©G(«-l)0Gm,m(,* J ; Where 

(1g ^ wg ) 

U^,M = (7^^^"^^""¥Ethg'®"'"^''""^|e2eG-)(e°l a^d 

UmA^) = (7r^'"''""¥Eti7!i'®"'"^°""^|eleGr.)(4l;^ : l-^V- Here, the sym- 
bol |e^) represents the fcth state of the ith MUB {i : 0, ...d). It is worth noting that 

, (eG((»-l)0G'"eG")0G'")U 

in odd prime power dimensions the phase factor ^ ,,° ,,„ — t: — ^— [ — ^ , ' — r— r is 

^ ^ ^ ( 0G(('-l)0G'"0G"'))2(-y®<3*"'®G"')2 

always equal to 1 so that the Wigner operators are invariant in all MUB's (up to a 
relabelling) j^OJ. The present treatment is also valid in even prime power dimensions. 

Let us now evaluate the value of the averaged sum of the Wigner operators 
(expressed in the computational basis) along a vertical straight line of the phase 
space: 

J-lsr^d-l TT/0 _ J-2x^d-l x^d-1 / 0G(n0GnffiG'n0GJ2)\ jj ( n.\ 

" Z^i2=0 '^'^(ii.Ja) ~^ 2^i2=0 2^m,n=0 UC ) ^m.,n{^) 

= d-'.Et-Lo rf.5n.,o(7^"^"'^""^) U^AO) = d-\Et-Lo (7^^^"®^"^) t/o,n(0) = 

V^rf-1 ^-1 /^0G(n0G'l)N/. ©G(fc0G")N |„0\/„0| _ W-1 Y^rf-1 ^S UO\/cO| 

Z^m,n,fc=0 " -WG J W G ) Pfc/\^fcl " " • Z^m,n,fc=0 "-"n,*: Pfc/V^fcl 

Let us finally evaluate the value of the averaged sum of the Wigner operators 
along a non-vertical straight line of the phase space (of slope k — l,k : l...d)] this 
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sum is equal to c?"^- Ea,m,n=o...d-i7G''^'®''"®''*'®''"'f^m,n(0) where ,i^ = ao + a,i2 = 
(fc-l)©Ga; 

we can rewrite it in the form 

^-2 Y^ eGi[@Gn'®Gl^@Gm' ^^(eG((''-l)@G^eG^)@G^))^ 

" •Z^/3,m,n=o...d-i7G • ^^eG((fc-i)©G™©G™))^ (^®g(™©g"))^ •'^m'.n'l'Cj wuere 

^l = (fc - 1) ©G ii 0G 12 = {k - 1) ©G «o, "^2 = ii = P,rn' = {k - 1) QcmQcn and 
ra' = -m. summing firstly over j3 and making use of the identity (|14p. we obtain a 

, (eG(('=-i)©G™eG")©G™)p 
factor (i.(5m/ n; now, when m' = then „ ,,^ ,,^ — p^ — r— [ — . . 1, — r— r = 1: the sum 

'"''-" ' (- Og(('=-i)©g'"©g'")\^C^®g(™©g"')\^ ' 

reduces thus to 

d'^- J2n'=o...d-i Ig^ ^ ' ^ -Uoyik); this is the sum of the Wigner operators along 
a vertical line, with the operators rewritten in the fcth MUB; in virtue of the previous 
result, it is equal to |ef^)(e*^J, the projector onto the iith state of the kth MUB. 
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Table 1: The field (Galois) multiplication in dimension 4. 
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Table 2: The field (Galois) addition in dimension 4. 



Appendix 2: Galois operations in dimension 4 
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